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Abstract 

Recent ultrafast experiments have implicated intrachain base-stacking rather than 
base-pairing as the crucial factor in determining the fate and transport of pho- 
toexcited species in DNA chains. An important issue that has emerged concerns 
whether or not a Frenkel excitons is sufficient one needs charge-transfer states to 
fully account for the dynamics. In we present an SU(2) <g> SU(2) lattice model 
which incorporates both intrachain and interchain electronic interactions to study 
the quantum mechanical evolution of an initial excitonic state placed on either the 
adenosine or thymidine side of a model B DNA poly (dA). poly (dT) duplex. Our 
calculations indicate that over several hundred femtoseconds, the adenosine exciton 
remains a cohesive excitonic wave packet on the adenosine side of the chain where 
as the thymidine exciton rapidly decomposes into mobile electron/hole pairs along 
the thymidine side of the chain. In both cases, the very little transfer to the other 
chain is seen over the time-scale of our calculations. We attribute the difference 
in these dynamics to the roughly 4:1 ratio of hole vs. electron mobility along the 
thymidine chain. We also show that this difference is is robust even when structural 
fluctuations are introduced in the form of static off-diagonal disorder. 

Key words: DNA, molecular biophysics, excited states, lattice theory, quantum 
chromodynamics 

PACS: 87.14.Gg, 87.15.Mi, 87.15.He, 31.15.Ar 



1 Introduction 



For all life-forms on Earth with the exception of certain viruses, genetic infor- 
mation is carried within the cellular nucleus via strands of strands of deoxyri- 
bonucleic acid (DNA). The genetic information itself is encoded in the specific 
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sequence of the nucleic acid bases: adenine (A), thymine (T), guanine (G), and 
cytidine (C). DNA is also a strong absorber of ultraviolet light leaving it highly 
susceptible to photomutagenic damage with the primary photoproducts be- 
ing bipyrimidine dimers linking neighboring T bases. For all organisms, this 
susceptibility is compensated for in part through enzymatic repair actions 
that remove damaged segments along one strand using the complementary 
strand as a template for replacement. Such repair mechanisms are quite costly 
energetically. Remarkably, however, DNA is intrinsically photochemically sta- 
ble as single bases are able to rapidly convert photoexcitation energy into 
thermal energy on a picosecond time scale through non-radiative electronic 
processes. What remains poorly understood, is the role that base-pairing and 
base-stacking plays in the transport and migration of the initial excitation 
along the double helix. Clearly, such factors are important since the UV ab- 
sorption of DNA largely represents the weighted sum of the absorption spectra 
of it constituent bases whereas the distribution of lesions formed as the result 
of photoexcitation are generally not uniformly distributed along the chain it- 
self and depend strongly upon sequence, suggesting some degree of coupling 
between bases. [1] 

Given the importance of DNA in biological system and its emerging role as a 
scaffold and conduit for electronic transport in molecular electronic devices, 
[2] DNA in its many forms is a well studied and well characterized system. 
What remains poorly understood, however, is the role that base-pairing and 
base-stacking plays in the transport and migration of the initial excitation 
along the double helix. [3,1] Such factors are important since the UV absorp- 
tion of DNA largely represents the weighted sum of the absorption spectra of 
it constituent bases whereas the distribution of lesions formed as the result 
of photoexcitation are generally not uniformly distributed along the chain it- 
self and depend strongly upon sequence, suggesting some degree of coupling 
between bases. [1] 

Recent work by various groups has underscored the different roles that base- 
stacking and base-pairing play in mediating the fate of an electronic excitation 
in DNA. [1,3] Over 40 years ago, Lowdin discussed proton tunneling between 
bases as a excited state deactivation mechanism in DNA[4] and evidence of this 
was recently reported by Schultz et al. [5] In contrast, however, ultrafast fluo- 
rescence of double helix poly (dA). poly (dT) oligomers by Crespo- Hernandez et 
al.[3] and by Markovitsi et al. [1] give compelling evidence that base-stacking 
rather than base-pairing largely determines the fate of an excited state in 
DNA chains composed of A and T bases with long-lived intrastrand states 
forming when ever A is stacked with itself or with T. However, there is con- 
siderable debate regarding whether or not the dynamics can be explained via 
purely Frenkel exciton models [6,7,8] or whether charge-transfer states play 
an intermediate role. [9] 
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Fig. 1. Three dimensional structure of stacked A-T base pairs along with the corre- 
sponding lattice model. 

Here we report on a series of quantum dynamical calculations that explore 
the fate of a localized exciton placed on either the A side or T side of the 
B DNA duplex poly(dA) 10 .poly(dT) 10 . Our theoretical model is based upon 
a SU(2) (g> SU(2) lattice model we recently introduced [10] that consists of 
localized hopping interactions for electrons and holes between adjacent base 
pairs along each strand (t a j) as well as cross-strand terms linking paired bases 
(hi) and "diagonal" terms which account for the 7r stacking interaction between 
base j on one chain and base j ± 1 on the other chain (rf) in which r~ 
denotes coupling in the 5'-5' direction and coupling in the 3'-3' direction. 
Fig. 1 shows the three-dimensional structure of poly(dA)i .poly(dT)i and 
the topology of the equivalent lattice model. We also consider here the role of 
geometric or structural fluctuations in the electronic dynamics. 



2 Theoretical Model 



Taking each link as Fig. 1 as a specific electron, hole, or excitonic, hopping 
term, we arrive at the following single particle Hamiltonian, 



j 

+^]+i(r+7+ + r"7_)^ + $(r+7+ + rj^^j+x, (1) 

where ipj and ipj are SU (2) spinors that act on the ground-state to create and 
remove an electron (or hole) on the jth adenosine or thymidine base along 
the chain. The 7 operators are the 2x2 Pauli spin matrices with ip ■ = 71^] 
and 7+ + 7_ = 71 providing the mixing between the two chains. As discussed 
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in Ref. [10], we used the highest occupied and lowest unoccupied (ir and n*) 
orbitals localized on each base as a orthonormal basis. For the single particle 
terms (representing electron and hole transfer between bases), we use values 
reported by Anatram and Mehrez as determined by computing the Coulomb 
integrals between HOMO and LUMO levels on adjacent base pairs with in 
a double-strand B DNA sequence using density functional theory (B3LYP/6- 
31G) [11] taking the geometries of each base from the B-DNA structure. When 
r+ = rj , Eq. 2 is identical to the Hamiltonian used by Creutz and Horvath [12] 
to describe chiral symmetry in quantum chromodynamics in which the terms 
proportional to r are introduced to make the "doublers" at q oc it heavier than 
the states at q oc 0. 

Taking the chain to homogeneous and infinite in extent, one can easily deter- 
mine the energy spectrum of the valence and conduction bands by diagonal- 
izing 



where e a ^ and t a ,b are the valence band or conduction band site energies and 
intra-strand hopping integrals. When rj~ = r~, Eq. 2 is identical to the Hamil- 
tonian used by Creutz and Horvath [12] to describe chiral symmetry in quan- 
tum chromodynamics in which the terms proportional to r are introduced to 
make the "doublers" at q oc 7r heavier than the states at q oc 0. In particular, 
we note that when t = h/2r, the band closes at q = ±n but has a gap at 



The single particle parameters are taken from Anatram and Mehrez as de- 
termined by computing the Coulomb integrals between HOMO and LUMO 
levels on adjacent base pairs with in a double-strand B DNA sequence using 
density functional theory (B3LYP/6-31G) [11]. Parameters used in our model 
are presented in Table 1. It is important to note thatg the asymmetry in- 
troduced with rlf 7^ rj gives directionality between the 3'- and 5'- ends of 
the chain. Introducing these parameters into Eq. 2 leads to 4 separate cosine 
shaped bands corresponding to conduction and valence bands localized along 
each chain as show in Fig. 2. 

The coupling between the conduction and valence bands is accomplished by in- 
troducing short-ranged Coulomb and exchange interactions as well as "dipole- 
dipole" terms which couple geminate electron-hole pairs on different sites and 
considering only single excitations, 



e a + 2t a cos(g) h + r + e~ iq + r~e +iq 
h + r + e +iq + r~e~ iq e b + 2t b cos(g) ) 



(2) 
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Fig. 2. (L) Band structure for homogeneous AT DNA chain. 



where the A m are spin-symmetrized composite operators that create or remove 
singlet or triplet electron/hole pairs in configuration |m) = \ij) where V mn = 
— (mn\\nm) + 25so(™^| \nm) where S = 1,0 is the total spin. [13,14,15,16] 
In our model, we include three types of electron/hole interactions, an on-site 
direct Coulomb J = (nn\\nn), an on-site exchange term K = (nn\nn), and an 
inter-site singlet exciton transfer term: d mn = 5so{mm\ \nn). 

For simplicity, we take the on-site Coulomb interaction as J = (nn\\nn) and 
the on-site exchange interaction K = (nn\nn) to be adjustable parameters 
with values J = 2.5 eV and K = 1.0 eV for both purines and pyrimidines. We 
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assume these interactions to be local since the distance at which the Coulomb 
energy between an electron/hole pair equals the thermal energy in aqueous 
ionic media at 300K is on the order of the base-stacking distance. Furthermore, 
such values are certainly in the correct range for conjugated cyclic organic sys- 
tems. [13,14,15,16] Lastly, we estimated the coupling between geminate elec- 
tron/hole pairs on different bases {nn\\mm) via a point-dipole approximation 
by mapping the tt — it* transition moments onto the corresponding base in 
the B DNA chain. These are obtained from the isolated bases by perform- 
ing single configuration interaction (CIS) calculations using the GAMESS[17] 
quantum chemistry package on the corresponding 9-methylated purines and 
1-methylated pyrimidines after optimizing the geometry at the HF/6-31(d)G 
level of theory. 

The use of the point-dipole approximation in this case is justified mostly 
for convenience and given the close proximity of the bases, multipole terms 
should be included in a more complete model. [18] As a result, the matrix 
elements used herein provide an upper limit (in magnitude) of the couplings 
between geminate electron-hole pairs. Most importantly, however, the point- 
dipole approximation provides a robust means of incorporating the geometric 
arrangement of the bases into our model. [10] 



3 Excited State Dynamics: Homogeneous Lattice 

We now consider the fate of an initial singlet electron/hole pair placed either 
in the middle of the thymidine side of the chain or the adenosine side of the 
chain (i.e. a localized exciton with the electron and hole starting on the same 
site). We assume that such a configuration is the result of an photoexcitation 
at the appropriate photon energy (4.87 eV for the thymidine exciton and 5.21 
eV for the adenosine exciton respecitively) and based upon the observation 
that the UV absorption spectra largely represents the weighted sum of the 
UV spectra of the constituent bases. Since these are not stationary states, 
they evolve according to the time-dependent Schrodinger equation, which we 
integrate using the Tchebychev expansion of the time-evolution operator. [19] 

In the top two frames of Fig. 3 we show the transient probability for finding 
an exciton placed on the adenosine (left) or thymidine (right) chain at time 
t = in some other excitonic configuration along either the adenosine or 
thymidine chain at some time t later. In both cases, negligible exciton density 
is transferred between chains and the excitons rapidly become delocalized and 
scatter ballistically down the DNA chain. 

There are some striking differences, however, between the exciton dynamics 
in adenosine versus those in thymidine. First, in comparing the d\\ matrix 
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Table 1 

Charge-transfer and exciton transfer terms for AT DNA. (B-FORM) 



Description Value Reference 



Intrachain e transfer: A{ — Ai+± 




tea — 


+0.024eV 


Ref. 


[11] 


Intrachain e transfer: T — Tj+i 




teb = 


-0.023eV 


Ref. 


[11] 


Interchain e transfer: Ai — Tj 




h e = 


+0.063 eV 


Ref. 


[11] 


Interchain 3'-3' e transfer: Ai — Tj + i 




rf = 


-0.012 eV 


Ref. 


[11] 


Interchain 5'-5' e transfer: — Tj_i 




r 7 = 


-0.016 eV 


Ref. 


[11] 


Intrachain /i transfer: Tj — Tj+i 




thb = 


-0.098eV 


Ref. 


[11] 


Intrachain /i transfer: Ai — Ai+i 




tha = 


+0.021eV 


Ref. 


[11] 


Interchain h transfer: Ai — Ti 




hh = 


+0.002eV 


Ref. 


[11] 


Interchain 3'-3' h transfer: Ai — Tj+i 




4 = 


-0.007 eV 


Ref. 


[11] 


Interchain 5'-5' h transfer: Ai — Tj_i 




r i = 


+0.050 eV 


Ref. 


[11] 


Intrachain parallel dipole-dipole: Tj — 


Ti+i 


4 = 


= 0.143eV 


Ref. 


[10] 


Intrachain parallel dipole-dipole: Ai - 


A i+ i 


df = 


0.0698eV 


Ref. 


[10] 


Interchain perpendicular dipole-dipole: Ai — T 


i d± ■- 


= -0.099 


Ref. 


[10] 


Interchain 3'-3' dipole-dipole: Ai — Ti 


fi 


d+ = 


-0.013eV 


Ref. 


[10] 


Interchain 5'-5' dipole-dipole: Ai — Ti 


-l 


dr = 


-0.006eV 


Ref. 


[10] 


Site energies 












A(LUMO) 






= 0.259eV 


Ref. 


[11] 


A(HOMO) 




£Ah = 


= -5.45eV 


Ref. 


[11] 


T(LUMO) 




ere = 


-0.931eV 


Ref. 


[11] 


T(HOMO) 




^Th = 


-6.298eV 


Ref. 


[11] 



elements, one easily concludes that the exciton mobility along the thymidine 
chain is considerably greater than the mobility along the adenosine chain. This 
can be see in Fig 2. comparing the time required for an excitonic wavepacket 
to reach the end of either chain. In adenosine, the exciton travels nearly 5 base 
pairs in about 25 fs where as an exciton along the thymidine chain covers the 
same distance in about 10 fs. This factor of two difference in the exciton ve- 
locity is commensurate with the rs1:2 ratio of the dj_ : d\ intrachain excitonic 
couplings. 

Secondly, we note that the adenosine exciton remains qualitatively more "co- 
hesive" than the thymidine exciton showing a number of ballistic traverses up 
and down the adenosine chain over the 300 fs we performed the calculation. 
One can also note that the exciton velocity in the 5'-3' direction is slightly 
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Fig. 3. Time evolution of excitonic density following excitation placed on T5 
(top-left) or A§ (top-right) site. The bottom panels show the net charge on each 
site following excitation of T5 (bottom-left) or A§ (bottom-right) 

greater than in the 3'-5' direction as evidenced by the exciton rebounds off 
site A\ slightly sooner than it rebounds from site A w . This is due to the asym- 
metry introduced in by the and terms. All in all, one can clearly note 
a series of strong recurrences for finding the adenosine exciton on the original 
site A§ every 100 fs. The thymidine exciton dynamics are far more complex 
as the exciton rapidly breaks apart. While few recursions can be noted, how- 
ever, after the first ballistic traverse, the thymidine exciton no longer exists 
as a cohesive wavepacket and is more or less uniformly distributed along the 
thymidine side of the chain. 

The excitonic dynamics only tell part of story. In the lower frames of Fig. 3 
we show the net charge taken as the difference between the hole density and 
electron density on a given base. In the case where the initial exciton is on the 
adenosine chain, very little charge-separation occurs over the time scale of our 
calculation. On the other hand, when the exciton is placed on the thymidine 



8 



Thymidine 



Adenosine 




0.8 
0.6 
0.4 
0.2 



20 40 60 80 100 



t(fs) 



20 40 60 80 100 



t(fs) 



Fig. 4. Net probability for an exciton originating on the thymidine (left) or adenosine 
(right) side of the chain to remain as an exciton on the original chain, to be found 
as an exciton on the other chain, and total exciton population. 



chain, the exciton almost immediately evolves into a linear combination of 
excitonic and charge-separated configurations. What is also striking is that 
in neither case do either the electron or hole transfer over to the other chain 
even though energetically charge-separated states with the electron on the 
thymidine and the hole on the adenosine sides of the chain are the lowest 
energy states of our model. [10] It is possible, that by including dissipation or 
decoherence into our dynamics, such relaxation will occur, however, on a time 
scale dictated by cross-chain transfer terms. For the coupling terms at hand, 
electron or hole transfer across base pairs occurs on the time scale of 3-4 ps. 

The difference between the excitonic dynamics following excitation of A vs. T 
can be quantitatively noted by comparing the curves shown in Fig. 4 where we 
compare the projection of the time-evolved state onto the excitonic configura- 
tions of the chain on which the exciton was placed (Paa an d Ptt) compared 
to the projection onto the excitonic configurations of the other chain {Pat and 
Pta)- For the case in which the adenosine chain was excited, approximately 
75% of the total probability density remains as excitonic configurations along 
the adenosine chain. In stark contrast, only about 40% of the initial thymidine 
exciton remains excitonic along the thymidine side of the chain. The reason 
for the remarkable difference between the two chains stems from the difference 
in electron and hole mobility along the thymidine chain. Indeed, comparing 
the electron and hole hopping terms given above, th/t e ~ 4 for the thymidine 
chain compared to t h /t e pa 1 for along the adenosine chain. This is manifest 
in the lower left panel of Fig. 3 where we see almost immediately a nega- 
tive charge remaining for a few fs on the site where the initial excitation was 
placed. In contrast, no charge-separation is seen following excitation along the 
adenosine side of the chain. 
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Fig. 5. Same as in Fig. 4 except averaging over ensemble of configurations. 
4 Excited State Dynamics: Effects of Static Disorder 

A proper model of the electronic dynamics in DNA must include some con- 
tribution from the solvent and environment. For DNA in water at 300K, we 
assume that the electronic processes described in our model are fast compared 
to the structural and environmental fluctuations the DNA lattice itself such 
that the parameters for the B-DNA structure represent the average values for 
the system in aqueous media. For example, the characteristic time scale for 
the relative lateral and longitudinal motions of bases in DNA is 10 to 100 fs 
with amplitudes of 0.01 to 0.1 nm. [20] Since electronic interactions between 
bases are sensitive to the fluctuations in the geometry of the DNA structure, 
factors such as salt concentration and other solvent media will have profound 
impact on the structure and hence on the model parameters. 

Since short strands of DNA are fairly rigid, the electronic coupling terms are 
likely most sensitive to the upon the dihedral angle, % between adjacent 
bases. If we take the fluctuations in 6^ to be 59 2 = kBT/IQ 2 where / is 
the reduced moment of inertia of the AT base-pair and Q = 25cm~ 1 is the 
torsional frequency. This gives an RMS angular fluctuation of about 5% about 
the avg. = 35.4° helical angle. Since this is a small angular deviation, 

we take the fluctuations in the electronic terms to be proportional to S9 2 and 
sample these terms from normal distributions about B-DNA average values. 

The effect of disorder on the excitonic dynamics are seen in Fig. 5. As in case 
of the homogeneous lattice above, our initial state is a Frenkel exciton on 
thymidine #5. The right-hand plot of Fig. 5 shows the probability of finding 
the time-evolved state in various excitonic configurations. While the time- 
scale for the exciton break up is slower than in the homogenous lattice, one 
can deduce that even with static off-diagonal disorder, the thymidine exciton 
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Fig. 6. RMS width of electron and hole densities vs. time for excitations originating 
on different sides of the AT chain. Top: Homogeneous chain, Bottom: averaging over 
20 random lattice configurations. 

is subject to dissociate to charge-transfer (i.e. polaron) pairs on an ultrafast 
timescale. The adenosine exciton appears to be more profoundly affected by 
the lattice disorder with some population transfered to the thymidine chain. 
Since the total exciton population is lowered upon introducing disorder, the 
remaining excited state population exists as charge-separated pairs. Analysis 
of the transient excited state wavepackets (not-shown) indicates a small but 
significant fraction of Af — Tf interstrand charge-transfer configurations. 

Lastly, we examine the relative mobilities of the electron and hole by calculat- 
ing the r.m.s. width of the electron and hole densities on either chain following 
excitation, the results of which are shown in Fig. 6. In the top two plots are the 
r.m.s. widths for the homogeneous lattice following excitation. As discussed 
above, the hole is considerably more mobile than the electron along the thymi- 
dine side of the chain. This is very much evidenced by the more rapid spread 
of the hole density compared to the the electron density. Likewise, along the 
adenosine side, the conduction and valence bands have nearly identical band- 
widths and the electron and hole densities evolve very much in concert. This 
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scenario is preserved upon even averaging over configurations indicating that 
to within a 5% error in our parameter set or the to within the typical thermal 
fluctuations of a DNA chain, our predicted dynamics are quite robust. 



5 Summary 

The results described herein paint a similar picture to that described by recent 
ultrafast spectroscopic investigations of (dA).(dT) oligomers in that the ini- 
tial excitonic dynamics is dominated by base-stacking type interactions rather 
than by inter-base couplings. Interchain transfer is multiple orders of magni- 
tude slower than the intrachain transport of both geminate electron/hole pairs 
as excitons and independent charge-separated species. Indeed, for an exciton 
placed on the adenosine chain, our model predicts that exciton remains as a 
largely cohesive and geminate electron/hole pair wave function as it scatters 
along the adenosine side of the chain. Our model also highlights how the dif- 
ference between the mobilities in the conduction and valence bands localized 
along each chain impact the excitonic dynamics by facilitating the break up 
of the thymidine exciton into separate mobile charge-carriers. In the actual 
physical system, the mobility of the free electron and hole along the chain will 
certainly be dressed by the polarization of the medium and reorganization of 
the lattice such that the coherent transport depicted here will be replaced by 
incoherent hopping between bases. 

The significance of the breakup of the exciton is twofold. First, it is well rec- 
ognized that photoexcitation of adjacent stacked pyrimidine bases leads to 
the formation of cis-syn cyclobutane pyrimidine-dimer lesions. However, this 
dimerization occurs only in the triplet (rather than singlet) excited state. [21] 
Consequently, spin-flip must occur ether via spin-orbit coupling or via recom- 
bination of polaron pairs. [16] If we assume that the spins are decorrelated at 
some intermediate distance r oc e 2 /ekT where the Coulomb energy is equal to 
the thermal energy, photoexcitation of a thymine sequence could rapidly result 
in a population of triplet excitons formed by exciton dissociation followed by 
geminate recombination. Secondly, the process is reversible and triplet reacti- 
vation of the dimer can lead to repair of the lesion. 

Isolating the photoexcitation to the originally excited chain minimizes the po- 
tential mutagenenic damage to the DNA sequence since it preserves the com- 
plementary chain as an undamaged back-up copy of the genetic information. 
It is fascinating to speculate whether or not the isolation of a photoexcita- 
tion and its photoproducts to the original chain was an early evolutionary 
selection criteria for the eventual emergence of DNA as the carrier of genetic 
information. 
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In conclusion, we present herein a rather compelling model for the short- 
time dynamics of the excited states in DNA chains that incorporates both 
charge-transfer and excitonic transfer. It is certainly not a complete model and 
parametric refinements are warranted before quantitative predictions can be 
established. For certain, there are various potentially important contributions 
we have left out: disorder in the system, the fluctuations and vibrations of 
the lattice, polarization of the media, dissipation, quantum decoherence. We 
hope that this work serves as a starting point for including these physical 
interactions into a more comprehensive description of this system. 
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